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ABSTRACT 


This  report  describes  a  base  technique  for  the  24-  and  48-hr 
prediction  of  stratospheric  contour  height  changes  in  winter  at  100,  50, 
and  30  mb.  On  independent  data,  this  technique  yields  superior  results  to 
persistence  at  all  three  levels  and  for  both  forecast  intervals.  Prediction 
equations  are  derived  by  applying  the  screening  regression  technique  to 
atmospheric  variables  at  a  network  of  grid  points  surrounding  a  predicand 
point.  Incorporation  of  predictors,  based  on  perfect  prognoses  at  lower  levels, 
brings  about  a  significant  improvement  in  the  results.  Some  improvement 
is  also  noted  when  a  geographical  stratification  is  employed.  However, 
orientation  of  the  grid  network  with  the  flow  pattern  did  not  result  in  any  substantial 
improvement. 
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SECTION  I 


INTRODUCTION 

The  objective  of  this  study  is  to  demonstrate  the  feasibility  of  stratospheric 
prediction  by  developing  a  physical-statistical  base  technique  to  forecast  100-,  50-, 
and  30-mb  heights  for  24-  and  48-hr  periods. 

Previous  work  in  this  area  was  aimed  at  developing  regression  equations  based 
on  middle-  and  upper -tropospheric  prognoses  to  extrapolate  upward  to  derive  temper¬ 
atures,  winds,  and  heights  at  stratospheric  levels  [6],  The  results  of  this  earlier  test 
were  not  positive  and  demonstrated  the  need  for  a  technique  superior  to  persistence. 

It  was  suggested  at  the  time  that  a  technique  using  a  more  direct  approacli  instead  of 
the  vertical-extrapolation  equations  should  be  investigated. 

The  framework  of  the  present  study  called  for  a  rather  modest  effort.  The 
experimental  design  required  a  minimum  amount  of  sophistication  and  a  restricted 
geographical  application.  The  so-called  “persistence”  technique  was  used  as  a 
control  to  demonstrate  the  feasibility  of  more  extensive  and  elaborate  investigations. 
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SECTION  II 


DATA  PROCESSING 

Hemispheric  grid-point  data,  provided  by  the  Air  Force  Global  Weather 
Central  (GWC)  and  containing  heights  at  six  constant  pressure  levels  determined 
twice  daily  for  December  1963  and  January  and  February  1964,  were  available  ior 
this  study. 

The  area  chosen  for  feasibility  testing  is  shown  in  Fig.  1.  The  48  predictand 
points  within  the  area  are  located  at  every  other  Joint  Numerical  Weather  Prediction 
(JNWP)  grid  point. 

Selection  of  cases  was  limited  to  the  81  consecutive  map  times  extending 
from  0000  GMT,  1  December  1963  to  0000  GMT,  12  January  1964,  yielding  3888 
cases  (48  cases  per  map  time).  Of  these,  16  maps  (768  cases)  were  made  avail¬ 
able  for  independent -data  verification  by  withholding  the  data  of  every  fifth  map 
time,  leaving  a  sample  of  3120  cases  for  development  work. 
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Fig.  1.  Area  used  for  feasibility  tests.  Grid  points  shown  are  predictand  points. 


SECTION  III 


THE  PREDICTION  TECHNIQUE 

1.  The  Composite  Grid 

A  grid  for  extracting  predictor  information  surrounding  any  predictand  point  was 
constructed  so  that  the  grid  would  measure  variables  at  locations  relative  to  a  predict¬ 
and  point  rather  than  at  fixed  geographical  locations  (See  Fig.  2).  The  grid  point  defined 
by  the  (K,L)-location  (3,3)  is  placed  over  the  predictand  point,  and  the  grid  is  oriented 
so  that  the  line  K=3  coincides  with  the  meridian  passing  through  the  predictand  point. 
For  development  work,  grid  placement  and  data  tabulation  were  done  by  computer  pro¬ 
grams,  and  “analyzed  maps’’  were  on  magnetic  tape  (an  option  in  the  computer  program 
permits  the  employment  of  an  alternative  grid  orientation  —  that  is,  with  respect  to  the 
100-mb  flow  rather  than  to  north-south).  On  a  polar  stereographic  projection  with 
standard  parallel  at  60°N,  the  5  X  5  array  forms  a  set  of  evenly-spaced  points  with 
the  grid  interval  being  equivalent  to  two  JNWP  grid  intervals  (762  km  at  60°N).  The 
25  points  defined  by  this  grid  system  were  the  ones  used  for  basic  predictor  tabulation. 

2.  Screening  Regression 

The  screening  procedure  suggested  by  Bryan  l  1]  and  developed  for  the  IBM  704 
electronic  computer  by  Miller  [3,  4]  was  used  to  screen  the  possible  predictors  identi¬ 
fied  in  subsequent  sections  (this  program  has  also  been  written  for  the  IBM  7094). 

One  who  designs  a  statistical  prediction  experiment  invariably  likes  to  consider  all 
predictors  deemed  important  on  the  basis  of  previous  theoretical,  synoptic,  and  em- 
pirical  work,  but  as  Lorenz  [2]  points  out,  a  prediction  equation  should  contain  few 
predictors  in  comparison  with  the  size  of  the  developmental  sample;  if  there  are  too 
many,  a  relationship  that  fits  the  sample  used  to  establish  it  is  likely  to  fail  when 
applied  to  a  new  sample.  The  object  of  the  screening  procedure  is  to  select  from  a 
set  of  possible  predictors  the  subset  that  most  significantly  and  independently  con¬ 
tributes  to  reducing  the  variance  of  the  predictand. 

From  an  array  of  possible  predictors,  the  screening  procedure  first  selects 
the  one  that  has  the  highest  linear  correlation  with  the  predictand  in  question.  This 
predictor  is  then  held  constant  and  partial-correlation  coefficients  between  the  pre¬ 
dictand  and  each  of  the  remaining  predictors  are  examined;  the  predictor  now  asso¬ 
ciated  with  the  highest  coefficient  is  the  second  one  selected.  Additional  predictors 
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Fig.  2.  Composite  grid  overlay.  Predictand  point 
is  at  K=3,  L=3.  Predictor  points  are  circled. 
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are  chosen  similarly.  Selection  is  halted  whenever  a  predictor  fails  to  pass  a  signifi¬ 
cance  test.  After  the  significant  predictors  have  been  selected,  the  regression  coeffi¬ 
cients  are  obtained  by  the  method  of  least  squares. 

The  criterion  of  significance,  as  applied  to  the  screening  procedure,  is  not  clear- 
cut  because  the  usual  F-test  methods  (e.g.,  [5])  are  not  applicable  [4],  If  a  predictor  is 
chosen  at  random  from  a  group  of  predictors,  an  F-test  is  usually  taken  at  the  95% 
level;  this  allows  a  l-in-20  chance  of  considering  the  predictor  significant  when,  in  fact, 
it  is  not.  Because  the  screening  procedure  does  not  select  predictors  randomly,  a  more 
severe  test  is  needed  to  specify  a  l-in-20  chance.  For  his  screening  procedure,  Miller 
[4]  suggested  that  the  critical  F-value  be  a  function  of  the  number  of  possible  predictors. 
The  F-test  was  used  in  this  form  in  these  experiments. 

3.  Predictands 

Height  changes  at  100,  50,  and  30  mb  for  forecast  intervals  of  24  and  48  hr  were 
chosen  as  predictands  for  the  study.  The  predictand  list  is  shown  in  Table  I. 

4.  Predictors  Considered 

The  GWC  hemispheric  grid-point  data  were  used  as  the  primary  source  of  predictor 
data.  Special  preprocessing  programs  automatically  derived  5X5  grid-point  arrays  of 
height  and  thickness  data  for  each  predictand  point  in  the  developmental  sample.  In 
addition,  various  vorticity  terms  were  computed  at  the  predictand  point  by  conventional 
finite-difference  procedures.  The  list  of  570  possible  predictors  is  given  in  Table  II. 

Because  Miller’s  screening-regression  technique  (as  programmed  for  the  IBM 
7094)  has  an  upper  limit  of  180  predictors  which  can  be  examined  simultaneously,  the 
number  of  possible  predictors  had  to  be  reduced  subjectively  before  screening  regression 
was  applied. 
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TABLE  I 

LIST  OF  PREDICTANDS 

Symbol 

Description 

a2100  (24) 

24-hr  forecast  100-mb  height  change* 

AZ50  (24) 

24-hr  forecast  50-mb  height  change* 

AZ30  (24) 

24-hr  forecast  30-mb  height  change* 

a2ioo  <48> 

48-hr  forecast  100-mb  height  change* 

AZ50  (48) 

48-hr  forecast  50-mb  height  change* 

AZ30  (48) 

48-hr  forecast  30-mb  height  change* 

*Unit  of  measure  =  10  ft. 
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SECTION  IV 


PREDICTION  EXPERIMENTS 

A  series  of  experiments  was  formulated  to  examine  various  alternative  approach¬ 
es  to  stratospheric  prediction  within  the  limitations  of  the  feasibility  study.  The  initial 
experiment  consisted  of  orienting  our  grid  system  north-south,  applying  it  over  the  en¬ 
tire  predictand  area  (no  stratification),  and  using  only  predictors  from  Table  II  which 
did  not  incorporate  prognostic  information.  This  was  our  “base-technique’’  experiment. 
Subsequent  experiments  were  devised  to  determine  the  advantages  of  incorporating 
lower-level  prognostic-type  predictors,  orienting  the  grid  system  with  respect  to  the 
flow  pattern,  and  developing  a  simple  geographical  stratification  scheme.  With  the  ex¬ 
ception  of  the  incorporation  of  prognoses,  all  experiments  used  the  same  set  of  possible 
predictors  from  the  list  in  Table  II. 

Because  the  number  of  possible  predictors  exceeded  the  screening  program’s 
limit  of  180,  it  was  necessary  to  subjectively  reduce  the  predictor  list  before  screen¬ 
ing.  In  this  study,  the  reduction  was  accomplished  by  considering  seven  predictors  in 
each  subset  of  25.  These  seven  points  are  shown  on  the  grid  overlay  in  Fig.  2.  Other 
combinations  of  predictors  were  not  attempted  in  these  initial  experiments.  The  types 
of  experiments  performed  are  outlined  in  Table  III. 
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TABLE  II 

POSSIBLE  PREDICTORS 

Symbol 

Description 

Z500 

500-mb  height* 

Z300 

300-mb  height* 

Z200 

200-mb  height* 

zioo 

100-mb  height* 

Z50 

50-mb  height* 

N 

CO 

o 

30-mb  height* 

AZ... 

500 

12-hr  500-mb  height  change* 

300 

12-hr  300-mb  height  change* 

AZ 

200 

12-hr  200-mb  height  change* 

AZ 

100 

12-hr  100-mb  height  change* 

AZ,_n 

50 

12-hr  50-mb  height  change* 

AZoa 

30 

12-hr  30-mb  height  change* 

H500-300 

500-  to  300-mb  thickness* 

H200-100 

200-  to  100-mb  thickness* 

H100-30 

100-  to  30-mb  thickness* 

AH500-300 

12-hr  500-  to  300-mb  thickness  change* 

AH200-100 

12-hr  200-  to  100-mb  thickness  change* 

AH 100-30 

12-hr  100-  to  30-mb  thickness  change* 

♦Unit  of  measure  is  10  ft;  number  available  is  25. 
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TABLE  II  (cont’d) 


Symbol 

Description 

aZ500  (24) 

24-hr  forecast  500-mb  height  change* 

AZ200  (24) 

24-hr  forecast  200-mb  height  change* 

AZ500  (48) 

48-hr  forecast  500-mb  height  change* 

AZ200  (48> 

48-hr  forecast  200-mb  height  change* 

^500 

500-  mb  absolute  vorticityt 

??300 

300-mb  absolute  vorticityt 

^200 

200-mb  absolute  vorticityt 

^100 

100-mb  absolute  vorticityt 

V50 

50-mb  absolute  vorticityt 

V30 

30-mb  absolute  vorticityt 

ATJ500 

12-hr  500-mb  vorticity  changet 

ATJ300 

12-hr  300-mb  vorticity  changet 

At?200 

12-hr  200-mb  vorticity  changet 

Ai?ioo 

12-hr  100-mb  vorticity  changet 

AT?50 

12-hr  50-mb  vorticity  changet 

t> 

CO 

o 

12-hr  30-mb  vorticity  changet 

£T 

b  500-300 

500-  to  300-mb  thermal  vorticityt 

ST200-100 

200-  to  100-mb  thermal  vorticityt 

*Unit  of  measure  =  10  ft;  number  available  =  25. 
tUnit  of  measure  =  10“®  sec“l;  number  available 
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TABLE  II  (cont’d) 


Symbol 

Description 

£T 

s  100-30 

100-  to  30-mb  thermal  vorticityt 

A£T 

55  500-300 

12-hr  500-  to  300-mb  thermal  vorticity 
change t 

A£T 

s  200-100 

12-hr  200-  to  100-mb  thermal  vorticity 
change t 

A£T 

S  100-30 

12-hr  100-  to  30-mb  thermal  vorticity 
change t 

V 

Magnitude  of  100-mb  geostrophic  wind? 

V2 

Square  of  magnitude  of  100-mb  geostrophic 
wind  IT 

fUnit  of  measure  =  10“5  sec“l;  number  available  =  1. 
tUnit  of  measure  =  knots;  number  available  =  1. 

ITUnit  of  measure  =  (knots )%;  number  available  =  1. 
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TABLE  III 

STRATOSPHERIC  PREDICTION  EXPERIMENTS 


Exp. 

no. 

Grid  orientation 

Predictors 

Stratification 

1 

North-south 

No  prognoses 

Unstratified 

2 

North-south 

With  prognoses 

Unstratified 

3 

Flow 

No  prognoses 

Unstratified 

4 

Flow 

With  prognoses 

Unstratified 

5 

North- south 

No  prognoses 

Stratified 

6 

North-south 

With  prognoses 

Stratified 

12 


SECTION  V 


RESULTS 

In  all  of  the  experiments  described  in  this  section,  the  screening  regression 
technique  was  applied  to  3120  cases  to  develop  prediction  equations  for  24-  and  48-hr 
height  changes  at  100,  50,  and  30  mb.  These  equations  were  then  applied  to  the  768  cases 
which  comprised  the  independent  data  sample  (the  various  equations  can  be  found  in  the 
appendix) . 

5.  Experiment  1 

In  this  experiment,  the  grid  is  oriented  north-south.  There  is  no  stratification 
and  the  possible  predictors  listed  in  Table  II  (excluding  lower-level  prognostic  pre¬ 
dictors)  are  used  to  derive  regression  equations.  Table  IV  lists  the  predictors  in  the 
order  of  their  selection  by  the  screening  procedure,  and  the  percentage  of  the  total  vari¬ 
ance  explained  by  each.  The  predictor  symbols  are  defined  in  Table  II  and  the  accompa¬ 
nying  numbers  refer  to  the  (K,L)-predictor  locations  in  the  grid  system  shown  in 
Fig.  2. 

From  Table  IV  it  can  be  seen  that  the  first  predictor  selected  is  usually  the  12-hr 
height  change  ‘‘upstream”  for  the  particular  level  in  question.  This  grid  point  is  two  grid 
intervals  to  the  west  of  the  predictand  point.  The  only  exception  is  the  24-hr,  30-mb 
height  change  forecast  (AZ3q),  where  the  first  predictor  selected  is  the  12-hr,  50-mb 
height  change  located  at  the  predictand  point  (3,3). 

6.  Experiment  2 

This  experiment  considers  the  addition  of  500-  and  200-mb  height  prognoses  as 
possible  predictors.  The  reason  for  this  kind  of  test  is  that,  for  operational  purposes, 
one  may  have  available  a  good  set  of  mid-tropospheric  and  lower-stratospheric  dynamic 
prognoses  which  could  be  used  for  predictive  information.  In  these  particular  tests,  it 
must  be  kept  in  mind  that  we  have  used  actual  analyses  (perfect  prognoses);  one  must 
still  examine  the  problem  of  how  much  skill  is  lost  in  going  from  a  perfect  prognosis  to 
an  operational  prognosis. 
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TABLE  IV 

PREDICTORS  SELECTED  BY  SCREENING  REGRESSION  FOR  EXP.  1 
(north-south  orientation,  no  prognoses,  unstratified) 

(a)  24-hr  forecast  interval 


Order  of 

selection 

AZ100 

A250 

30 

predictor 

% 

reduction 

predictor 

% 

reduction 

predictor 

% 

reduction 

1 

AZ100(1-  3> 

23.4 

AZ50(1’  3) 

22.3 

AZ50<3’  3> 

21.6 

2 

AZ500<3’  3> 

5.6 

AZioo<3’  3> 

5.0 

AZ30(1-  3) 

5.9 

3 

AW3’3) 

1.7 

AH100-30<3’ 3> 

2.9 

H100-30<3’  3) 

2.9 

4 

200-100*  * 

5) 

1.6 

AZ30(1’  5) 

1.9 

H200-100(3'3) 

2.5 

5 

CT100-30(3, 

3) 

1.5 

H200-100(5,  5) 

1.7 

A”50(3'3) 

1.6 

6 

AZ30(1-  5) 

1.2 

A"50(3’3) 

1.5 

AT)30(3’3) 

1.3 

7 

H200-100(1, 

1) 

1.0 

AZ500(1’  5) 

0.9 

H100-30(1,  5) 

1.1 

8 

zioo<3'  3> 

1.2 

Z500(1’ 

0.8 

AZ50(1>  3) 

0.8 

9 

Z100(1’  3> 

1.8 

Z50<3’  3) 

1.7 

Z50(5’  5) 

0.7 

10 

Z30(5’  1J 

1.4 

Z30(5> 

0.8 

Z30(5’  1} 

2.4 

11 

AZ50(3’  3> 

1.4 

Z30(1’  5) 

1.7 

Z50(5> 

1.0 

12 

Z500<5’  3> 

1.0 

Z500(1'  5) 

0.9 

"50(3*  3) 

0.9 

13 

Z50(3’3) 

0.7 

H500-300(1>  5> 

0.7 

zioo<3’  3> 

1.7 

14 

— 

— 

AIJ30(3’  3) 

0.7 

— 

- 

15 

— 

— 

A£T100-30(3’  3) 

0.4 

— 

16 

- 

- 

AZ50(1- 

0.2 

— 

17 

— 

- 

AZ30(1’  3> 

0.3 

— 

— 

18 

- 

— 

H100-30<5’  3> 

0.3 

- 

— 

19 

— 

— 

z30(1.  3) 

0.5 

— 

— 

Total 

- 

43.5 

45.2 

- 

44.4 
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TABLE  IV  (cont’d) 


(b)  48-hr  forecast  interval 


Order  of 
selection 

AZ,„„ 

100 

50 

a2 

30 

predictor 

% 

reduction 

predictor 

% 

reduction 

predictor 

% 

reduction 

1 

AZioo(1’  3> 

17.5 

AZ50(1>  3) 

17.7 

AZ30(1’  3) 

19.9 

2 

AZ50(1’  5) 

3.1 

AZ30(1-  5) 

6.4 

H100-30<3-  3> 

8.5 

3 

H200-100(1'  5) 

2.8 

H100-30(1>  3> 

4.4 

AZ30(1’  5) 

3.1 

4 

zioo<3’  3) 

1.2 

z50(5,5)  . 

2.4 

AZ50(3'  3) 

2.5 

5 

O30(3,  3) 

7.7 

Z200(1'  X> 

2.0 

Z50(5’  5) 

1.9 

6 

H500-300(1’  5> 

3.3 

Z50<3'  3) 

1.9 

zioo(1- 

2.2 

7 

Z50(1’ 

2.5 

Z30<5’  A> 

1.5 

Z200<3’  3) 

2.4 

8 

H100-300(5,  ^ 

1.4 

Z50(1>  5) 

2.8 

Z30(5* 

1.3 

9 

Z500(5'  3> 

0.8 

Z500(1’  5> 

2.4 

H100-30(1-  5> 

2.4 

10 

AZ200(3’  3) 

1.4 

H200-100(1,  ^ 

2.1 

Z50<5’ 

1.4 

11 

W- 

1.1 

— 

- 

W®’  3) 

1.2 

12 

AZ50(1-  1} 

0.7 

— 

— 

H500-300(1,  5) 

0.7 

13 

AH200-100(1-  3) 

0.4 

— 

— 

Z300  (1’  5> 

2.1 

14 

Z300  (1'  5) 

0.4 

— 

— 

Z50  (1>  3> 

0.8 

15 

Z30  (1>  3) 

0.4 

— 

— 

^50  (3’  3) 

0.5 

16 

H100-30  (5’  5) 

0.4 

- 

— 

- 

— 

17 

Z300  (3’  3) 

0.4 

— 

— 

— 

— 

18 

Z50  (3’  3) 

0.4 

— 

- 

— 

— 

19 

Z200  5> 

0.4 

— 

- 

— 

— 

20 

H100-30  <3’  3> 

0.5 

— 

— 

— 

— 

Total 

- 

46.8 

— 

43.6 

~ 

50.9 
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The  first  predictors  selected  in  this  experiment  were  the  24-  and  48-hr,  200-mb 
prognostic  heights  at  the  predictand  point  for  the  24-  and  48-hr,  100-mb  height  predic¬ 
tion,  respectively  (see  Table  V).  At  50  mb,  the  first  predictors  selected  are  the  same 
as  in  Exp.  1,  with  the  prognostic  predictors  being  selected  second.  At  30  mb,  the  24-hr 
prediction  equation  does  not  select  a  prognostic  predictor  until  after  five  predictors, 
based  on  observed  data,  have  been  chosen,  while  for  the  48-hr  prediction  the  third  pre¬ 
dictor  selected  is  a  prognostic  one.  It  is  reasonable  to  expect  the  500-  and  200-mb 
prognostic  predictors  to  make  a  more  significant  contribution  at  the  lowest  (100-mb) 
level.  The  total  percent  reduction  of  variance  (PR)  for  this  experiment  is  higher  for 
each  of  the  six  predictands  than  the  corresponding  predictands  of  Exp.  1,  with  the 
difference  decreasing  with  increasing  height. 

7,  Experiment  3 

The  same  set  of  possible  predictors  used  in  Exp.  1  was  used  in  this  experiment. 
The  only  difference  was  in  the  selection  of  grid  orientation.  Whereas  predictors  in 
Exp.  1  were  derived  in  a  north-south  grid  orientation,  the  grid  for  this  experiment  was 
oriented  so  that  a  line  defined  by  K  =  3  (see  Fig.  2)  was  normal  to  the  100-mb  geostro- 
phic  wind  computed  at  the  predictand  point.  The  predictors  in  their  order  of  selection 
are  shown  in  Table  VI.  Note  that  the  first  predictor  selected  corresponds  to  that  of 
Exp.  1  (Table  IV)  for  all  six  predictands.  However,  although  the  coordinate  locations 
are  the  same,  the  geographical  locations  differ  because  of  the  difference  in  grid  orieta- 
tion.  While  the  percent  reduction  of  variance  attributed  to  these  flow-oriented  first 
predictors  is  greater  for  each  of  the  six  predictands,  the  total  PR  is  only  greater  for 
one  of  them  —  the  24-hr,  100-mb  height  change. 

8.  Experiment  4 

For  this  experiment,  the  grid  is  oriented  with  the  flow  and  the  predictor  list  is 
expanded  to  include  500-  and  200-mb  prognoses.  The  results  of  applying  the  screening 
procedure  are  shown  in  Table  VII.  Comparison  with  Exp.  2  (Table  V),  where  the  pre¬ 
dictor  list  is  the  same  but  the  orientation  is  different,  shows  that  here,  too,  one  selects 
the  same  predictor  first  for  all  six  predictands.  There  is  very  little  difference  between 
the  two  experiments  (2  and  4)  in  the  PR  for  the  first  selected  predictors,  the  flow- 
orientation  PR  was  generally  higher,  although  the  total  PR  for  all  six  predictands  was 
higher  for  the  north-south  orientation. 
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TABLE  V 

PREDICTORS  SELECTED  BY  SCREENING  REGRESSION  FOR  EXP.  2 
(north-south  orientation,  with  prognoses,  unstratified) 

(a)  24-hr  forecast  interval 


Order  of 
selection 

AZ 

100 

AZ_ 

50 

AZ 

30 

predictor 

% 

reduction 

predictor 

% 

reduction 

predictor 

% 

reduction 

1 

AZ200(24)  (3,3) 

56.6 

azc„ 

50 

(1,3) 

22.3 

AZCrt 

50 

(3,3) 

21.6 

2 

AZrn 

50 

(1,3) 

6.3 

AZ124) 

200 

(3,3) 

8.3 

AZ 

30 

(1,3) 

5.9 

3 

'*200-100 

(3,3) 

3.6 

AZoa 

30 

(3,3) 

5.2 

**100-30 

(3,3) 

2.9 

4 

**200-100 

(1,3) 

2.4 

**100-30 

(3,3) 

3.5 

**200-100 

(3,3) 

2.5 

5 

’*100-30 

(3,3) 

1.5 

AZ200(24) 

(5,3) 

2.6 

A"50 

(3,3) 

1.6 

6 

H100-30 

(1,3) 

1.7 

A**100-30 

(1,5) 

1.8 

AZ200<24)  (3>  3> 

1.5 

7 

AZr„ 

50 

(3,3) 

1.4 

An 

'50 

(3,3) 

1.3 

**100-30 

(1,5) 

2.0 

8 

AZ200(24)(5’3) 

0.6 

AZ„„ 

30 

(1,3) 

0.9 

AT)30 

(3,3) 

1.1 

9 

’*100-30 

(5,5) 

0.5 

AZ„f24) 
200v  ' 

(5,5) 

0.7 

Z50 

(5,5) 

1.0 

10 

A**100-30 

(1,5) 

0.5 

AZ  (24) 

500'  ’ 

(5,1) 

0.7 

Z30 

(5, 1) 

2.5 

11 

A^50 

(3,3) 

0.3 

**200-100 

(1,3) 

0.5 

az500(24)(5.D 

1.0 

12 

AZ200(24)(1’1) 

0.3 

Z500 

(1,1) 

0.6 

Z500 

(3,3) 

1.1 

13 

AZ500(24)  *5,  ■*> 

0.4 

Z50 

(3,3) 

1.2 

n300 

(3,3) 

0.8 

14 

AZ500 

(1,1) 

0.3 

**500-300 

(3,3) 

1.4 

Z50 

(5, 1) 

0.7 

15 

AZ500 

(3,3) 

0.3 

AZ500<24> 

(1,  1) 

1.2 

**200-100 

(1,3) 

0.8 

16 

**500-300 

(1,5) 

0.2 

Z50 

(1,5) 

0.7 

AZ500(24)  (*■’ 4) 

0.6 

17  . 

Z500 

(1,5) 

0.5 

Z30 

(5, 1) 

1.3 

Z500 

(1,1) 

0.8 

18 

Z500 

(3,3) 

0.3 

Z500 

(1,5) 

1.1 

AZ500(24)  <5’ 3) 

0.6 

19 

Z50 

(1,5) 

0.8 

A7J30 

(3,3) 

0.4 

AZioo 

(5,5) 

0.6 

20 

— 

- 

**200-100 

(1,5) 

0.4 

AZ500 

(5,3) 

0.4 

Total 

— 

78.5 

— 

56.1 

— 

50.0 
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TABLE  V  (cont’d) 


(b)  48 -hr  forecast  interval 


Order  of 
selection 

AZ 

100 

AZ 

50 

AZ_ 

30 

predictor 

% 

reduction 

predictor 

% 

reduction 

predictor 

% 

reduction 

1 

AZ200<48> 

(3,3) 

64.9 

AZr„ 

50 

(1,3) 

17.7 

AZ 

30 

(1,3) 

19.9 

2 

H200-100 

(3,3) 

6.7 

£Zn(±8) 
200'  ' 

(3,3) 

13.4 

H 100-30 

(3,3) 

8.5 

3 

AZ50 

(1.5) 

3.4 

H100-30 

(3,3) 

7.0 

AZ  (48) 

200'  ' 

|(3,3) 

5.9 

4 

AZ200(48) 

(5,3) 

2.7 

AZ„„ 

30 

(1,5) 

5.1 

AH100-30 

(1,5) 

3.8 

5 

AZr„ 

50 

(1,3) 

1.2 

AZ„  (48) 
200'  ' 

(5,3) 

4.3 

AZ50 

(3,3) 

3.2 

6 

H200-100 

(1,3) 

1.0 

AZ30 

(1,3) 

1.7 

AZ500(48) 

1  (5, 1) 

2.4 

7 

H 100-30 

(3,3) 

0.7 

AZ500<48> 

(5,1) 

1.3 

Z500 

(1,5) 

1.8 

8 

Z500 

(3,3) 

0.6 

AZ500(24) 

(1,1) 

1.3 

H500-300 

(1,5) 

1.8 

9 

Z300 

(1,3) 

0.8 

AZ200 

(1,3) 

1.0 

Z500 

(5,5) 

1.2 

10 

AZ500(48)(3’3) 

0.7 

Z500 

(1,5) 

0.8 

Z30 

(5,1) 

1.5 

11 

Z30 

(5,1) 

0.7 

Z300 

(1,5) 

1.8 

AZ500<48> 

(5,3) 

0.9 

12 

H200-100 

(1,5) 

0.4 

Z500 

(1,1) 

1.1 

H200-100 

(1,3) 

0.7 

13 

H500-300 

(1,5) 

0.5 

Z50 

(3,3) 

1.8 

AZ  (24) 

500'  ; 

(1,1) 

0.6 

14 

Z500 

(1,5) 

0.3 

Z50 

(5, 1) 

0.8 

H100-30 

(1,1) 

1.0 

15 

AZ500<48> 

(5,1) 

0.3 

Z50 

(1,5) 

1.7 

Z50 

(5,3) 

1.5 

16 

AZ500(24> 

(1,1) 

0.4 

H500-300 

(3,3) 

0.8 

H100-30 

(1,5) 

0.8 

17 

- 

— 

zioo 

(5,5) 

1.1 

Z500 

(3,3) 

1.1 

18 

- 

— 

zioo 

(5,3) 

0.9 

AZ50 

(1,5) 

0.7 

19 

— 

— 

H200-100 

(1,3) 

1.0 

AZ50 

(5,3) 

0.6 

20 

— 

- 

- 

- 

Z50 

(5,1) 

0.4 

Total 

- 

85.3 

- 

64.6 

- 

58.4 
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TABLE  VI 

PREDICTORS  SELECTED  BY  SCREENING  REGRESSION  FOR  EXP.  3 
(north-south  orientation,  no  prognoses,  unstratified) 

(a)  24-hr  forecast  interval 


Order  of 
selection 

AZ 

100 

a2 

50 

AZ 

30 

predictor 

% 

reduction 

predictor 

% 

reduction 

predictor 

% 

reduction 

1 

AZ100 

(1,3) 

26.2 

azc„ 

50 

(1,3) 

24.7 

AZ50 

(3,3) 

21.7 

2 

AZ500 

(3,3) 

5.3 

AZ 

100 

(3,3) 

4.8 

AZ30 

(1,3) 

9.2 

3 

Z30 

(1,3) 

2.6 

AZ 

30 

(1,3) 

3.0 

H100-30 

(3,3) 

2.7 

4 

AZ30 

(3,3) 

2.1 

AH100-30 

(3,3) 

1.6 

H200-100 

(3,3) 

2.6 

5 

Z30 

(5,3) 

2.0 

Z30 

(1,3) 

1.8 

Z30 

(5,1) 

1.3 

6 

Z50 

(3,3) 

1.5 

zioo 

(5,1) 

2.1 

A\o 

(3,3) 

1.2 

7 

Z500 

(5,5) 

1.4 

Ar?50 

(3,3) 

1.3 

A^30 

(3,3) 

0.9 

8 

A,,500 

(3,3) 

0.8 

Z50 

(3,3) 

1.1 

Z200 

(1,5) 

1.0 

9 

AZ50 

(1,3) 

0.7 

— 

— 

AZ„, 

30 

(5,3) 

0.5 

10 

zioo 

(3,3) 

1.0 

— 

— 

AZ50 

(1,3) 

0.5 

11 

Z300 

(5,3) 

3.4 

— 

— 

Z50 

(5,1) 

0.5 

12 

Z500 

(5,1) 

0.9 

— 

— 

H100-30 

(1,5) 

0.6 

13 

H500-300 

(3,3) 

0.5 

— 

- 

zioo 

(3,3) 

0.7 

14 

£T500-300(3,3^ 

0.6 

— 

— 

- 

— 

15 

AZ200 

(1,3) 

0.5 

— 

— 

— 

— 

16 

Z500 

(1,3) 

0.5 

— 

— 

— 

- 

17 

H500-300 

(1,5) 

0.7 

— 

— 

— 

— 

Total 

— 

50.7 

— 

40.4 

— 

43.4 
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TABLE  VI  (cont’d) 


(b)  48-hr  forecast  interval 


Order  of 
selection 

AZioo 

AZ50 

a2 

30 

predictor 

% 

reduction 

predictor 

% 

reduction 

predictor 

% 

reduction 

1 

AZ 

100 

(1,3) 

20.5 

AZ 

50 

(1,3) 

21.3 

AZn„ 

30 

(1,3) 

22.2 

2 

H200-100 

(1,5) 

5.0 

H100-30 

(1,3) 

5.4 

H100-30 

(3,3) 

8.3 

3 

AZ50 

(1,5) 

1.7 

AZ30 

(1,3) 

3.8 

AZrn 

50 

(3,3) 

3.1 

4 

zioo 

(3,3) 

1.5 

AZ30 

<1,  5) 

1.6 

H100-30 

(5,3) 

1.9 

5 

”30 

(3,3) 

6.1 

Z50 

(3,3) 

1.7 

AZ30 

(1,5) 

1.6 

6 

Z30 

(5,3) 

2.5 

zioo 

(5,1) 

3.4 

Z500 

(3,3) 

1.5 

7 

Z500 

(5,3) 

2.3 

Z30 

(5,3) 

0.9 

zioo 

(5,1) 

2.7 

8 

AZ50 

(1,3) 

1.1 

Z500 

(1,3) 

0.7 

A7,50 

(3,3) 

0.6 

9 

AH200-100(3’3) 

0.6 

Z500 

(5,3) 

1.0 

Z500 

(1,5) 

0.6 

10 

Z500 

(5,1) 

0.5 

AZ30 

(3,3) 

0.5 

'’soo 

(3,3) 

0.5 

11 

Z500 

(5,5) 

0.4 

Z50 

(1,3) 

0.5 

H200-100 

(1,1) 

0.6 

12 

AZ200 

(1,3) 

0.4 

AZr„ 

50 

(1,5) 

0.5 

AZ 

100 

(3,3) 

0.5 

13 

H500-300 

(3,3) 

0.5 

AZ500 

(1,3) 

0.2 

AZ 

50 

(1,3) 

0.5 

14 

Z200 

(1,5) 

0.4 

Z300 

(3,3) 

0.4 

H200-100 

(1,3) 

0.7 

15 

Z300 

(5,  1) 

0.3 

— 

- 

H500-300 

(5,1) 

0.3 

16 

Z200 

(5, 1) 

0.5 

— 

- 

Z200 

(3,3) 

0.4 

17 

zioo 

(5,3) 

0.5 

— 

- 

H500-300 

(1,3) 

0.4 

18 

- 

— 

- 

- 

Z300 

(3,3) 

0.4 

19 

— 

— 

— 

— 

zioo 

(5,3) 

0.4 

20 

— 

— 

— 

- 

Z50 

(5,3) 

0.5 

Total 

- 

44.8 

- 

41.9 

- 

47.7 
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TABLE  VII 

PREDICTORS  SELECTED  BY  SCREENING  REGRESSION  FOR  EXP.  4 
(flow  oriented,  with  prognoses,  unstratified) 

(a)  24-hr  forecast  interval 


Order  of 
selection 

AZ 

100 

AZe/1 

50 

AZon 

30 

predictor 

% 

reduction 

predictor 

% 

reduction 

predictor 

% 

reduction 

1 

AZ200(24)(3’3) 

56.5 

AZr„ 

50 

(1,3) 

24.7 

50 

(3,3) 

21.7 

2 

50 

(1,3) 

7.5 

f24) 
200v  1 

(3,3) 

6.5 

AZ_ 

30 

(1,3) 

9.2 

3 

ZiZ  (24) 
200' 

(5,3) 

3.5 

AZ_ 

30 

(3;  3) 

5.3 

H100-30 

(3,3) 

2.7 

4 

£T200-100(3,3) 

2.2 

H 100-30 

(3,3) 

3.6 

H200-100 

(3,3) 

2.6 

5 

H100-30 

(1,5) 

0.8 

AZ0  *24) 
200v  ' 

(5,3) 

2.6 

Z30 

(5, 1) 

1.3 

6 

AZ„„ 

50 

(3,3) 

0.8 

AZ30 

(1,3) 

1.8 

A7,50 

(3,3) 

1.2 

7 

AZioo 

(1,3) 

0.6 

A7?50 

(3,3) 

1.1 

A7,30 

(3,3) 

0.9 

8 

AH100-30 

(3,3) 

0.5 

_(24) 
200'  1 

(5,5) 

0.9 

Z200 

(1,5) 

1.0 

9 

H200-100 

(3,3) 

0.5 

AH200-100 

(1,3) 

0.7 

AZ200<24> 

(3,3) 

0.8 

10 

H 100-30 

(3,3) 

1.0 

— 

— 

H200-100 

(1,3) 

0.9 

11 

H200-100 

(1,3) 

0.7 

— 

- 

n300 

(3,3) 

0.7 

12 

H 100-30 

(1,3) 

0.9 

— 

- 

Z500 

(3,3) 

1.0 

13 

Z500 

(3,3) 

0.4 

— 

— 

H200-100 

(5,  5) 

0.6 

14 

Z300 

(5,3) 

0.8 

— 

- 

50 

(5,5) 

0.8 

15 

AZ500<24> 

(3,3) 

0.4 

— 

— 

— 

— 

16 

Z30 

(5,3) 

0.4 

- 

— 

— 

— 

17 

AZ200(24><1’1> 

0.3 

— 

— 

- 

— 

Total 

— 

77.8 

— 

47.2 

— 

45.4 
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TABLE  VII  (cont’d) 


(b)  48-hr  forecast  interval 
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9.  Experiments  5  and  6 


These  two  experiments  incorporated  the  idea  of  geographical  stratification  in  a 
simple,  straightforward  manner.  The  predictand  area  (Fig.  1)  was  merely  subdivided 
into  four  regions  (NW,  NE,  SE,  SW),  each  consisting  of  12  predictand  points  and  a  sample 
size  of  780  cases,  or  one-fourth  the  total  dependent  sample.  The  prognostic  predictors 
were  excluded  in  Exp.  5  and  included  in  Exp.  6  with  the  grid  oriented  north-south  for 
both  experiments.  The  predictors  in  the  order  of  their  selection  are  not  shown,  but  the 
total  PR’s  are  shown  in  Table  VIII,  which  summarizes  the  results  of  all  the  experiments 

on  dependent  data.  Generally  speaking,  the  first  predictor  selected  was  similar  to  those 
in  the  unstratified  experiments. 

10.  Independent  Data  Tests 

Equations  developed  from  the  six  types  of  experiments  described  above  were  applied 
to  a  set  oi  independent  data  consisting  of  768  cases  taken  from  16  map  times.  In  addition, 
the  technique  of  persistence  was  applied  to  the  sample  for  control  purposes  (the  state-of- 
the-ai  t  justiiies  this  kind  of  comparison  especially  at  50  and  30  mb).  Root-mean-square 
(ims)  errors  (in  feet)  are  shown  in  Table  IX.  For  Exps.  5  and  6,  the  results  for  the 
individual  four  stratified  areas  have  been  pooled  for  comparative  purposes.  Many  of  the 
comments  on  dependent  data  comparisons  apply  as  well  to  the  independent  data. 

The  introduction  of  prognostic  predictors  (Exp.  1  vs  2,  Exp.  3  vs  4,  Exp.  5  vs  6) 
results  in  significant  improvement,  particularly  at  100  mb.  This  is  due  to  the  rather 
high  correlation  between  height  changes  at  200  mb  (the  highest  predictor  level)  and  100  mb 
(the  lowest  predictand  level).  The  comparisons  between  north-south  and  flow  orientation 
(Exp.  1  vs  3,  Exp.  2  vs  4)  seem  to  indicate  that  there  is  little  to  be  gained  by  orienting 
the  grid  with  respect  to  the  flow.  The  crude  stratification  technique  employed  (Exp.  1  vs 
5,  Exp.  2  vs  6)  appears  to  have  been  successful,  especially  for  the  48-hr  forecast  interval. 
The  technique  of  persistence  yielded  the  largest  rms  errors  for  all  levels  and  forecast 
intervals. 

The  application  of  Exps.  1  and  2,  48-hr  prediction  equations,  and  persistence  to 
one  of  the  independent  data  situations  (1200  GMT,  30  December  1963)  is  shown  in 
Figs.  3-6.  The  superimposed  error  fields  (in  tens  of  feet)  are  represented  by  dashed 
lines.  Figure  3  shows  the  initial  100-mb  chart  of  1200  GMT,  30  December  1963.  The 
major  feature  of  this  map  is  the  large  trough  extending  from  east  of  Hudson  Bay  south- 
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RESULTS  ON  DEPENDENT  DATA 
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Fig.  3.  Observed  100  mb  height  analysis  for  1200  GMT,  30  December  1963.  Heights  are  in  hundreds  of 
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Fig.  4.  Experiment  1— 48-hr  100-mb  height  prognosis;  valid  1200  GMT,  1  January  1964.  Heights  are  in 
hundreds  of  feet.  Error  field  (in  tens  of  feet)  is  superimposed  in  dashed  lines. 
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Fig.  5.  Experiment  2— 48-hr  100-mb  height  prognosis;  valid  1200  GMT,  1  January  1964.  Heights  are  i 
hundreds  of  feet.  Error  field  (in  tens  of  feet)  is  superimposed  in  dashed  lines. 
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Fig.  6.  Observed  100-mb  height  analysis  for  1200  GMT,  1  January  1964.  Heights  are  in  hundreds  of 
feet.  Dashed  lines  are  the  error  field  (tens  of  feet)  for  48-hr  persistence. 


westward  to  Texas.  During  the  ensuing  48  hr,  the  northern  portion  of  this  trough 
moved  rapidly  eastward  leaving  a  rather  weak  trough  over  the  Mississippi  Valley 
area  (see  Fig.  6).  The  accompanying  48 -hr  height  changes  are  equal,  but  of  opposite 
sign,  to  the  error  field  of  the  persistence  technique  which  is  represented  in  Fig.  6  by 
the  dashed  lines.  Note  the  1400-ft  height  falls  in  the  upper  right  portion  of  the  map, 
while  most  of  the  remainder  of  the  map  is  characterized  by  height  rises,  with  a  height 
rise  center  in  excess  of  1200  ft  in  the  vicinity  of  James  Bay.  The  Exp.  1  prediction 
equation  results  (north-south  orientation,  no  prognoses,  unstratified)  are  shown  in 
Fig.  4.  There  has  been  a  broadening  in  the  trough  predicted  rather  than  an  eastward 
displacement.  The  main  problem  here  is  the  1400-ft  errors  near  the  upper  right  cor¬ 
ner  of  the  map.  Some  of  the  height  rises  over  the  continent  have  been  indicated  to 
some  extent  with  the  largest  errors  being  around  600  ft. 

The  results  using  the  Exp.  2  equation  (north-south  orientation,  with  prognoses) 
show  an  appreciable  improvement  (Fig.  5).  The  main  features  have  been  explained 
fairly  well:  the  eastward  displacement  of  the  northern  portion  of  the  trough,  the  linger¬ 
ing  trough  over  the  central  U.S.,  and  the  building  ridge  over  the  western  U.S.  The 
1200-ft  rises  over  James  Bay  have  been  almost  entirely  predicted  in  this  case.  The 
overall  rms  error  for  this  example  is  705  ft  for  persistence,  469  ft  for  Exp.  1,  and 
249  ft  for  Exp  2. 
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SECTION  VI 


CONCLUSIONS  AND  RECOMMENDATIONS 

The  successful  results  of  this  limited  feasibility  test  on  stratospheric- circulation 
prediction  clearly  indicate  that  further  research  of  greater  proportion  in  this  area  is 
warranted.  While  the  application  of  real  prognoses  may  not  yield  the  spectacular  re¬ 
sults  that  perfect  prognoses  do,  it  should  be  noted  that  even  the  simplest  base  technique 
devised  (north-south  orientation,  no  prognoses,  and  no  stratification)  yielded  results 
superior  to  persistence  for  all  three  levels  and  both  forecast  intervals.  Although  the 
use  of  a  flow  orientation  failed  to  yield  much  improvement,  the  adoption  of  a  simple 
stratification  scheme  gave  encouraging  results. 

A  logical  follow-on  research  plan  would  include: 

(a)  the  extension  of  technique  development  to  the  entire  Northern 
Hemisphere  and  to  other  seasons, 

(b)  derivation  of  additional  predictors;  for  example,  the  concept  of 
vorticity  conservation  could  be  employed  as  is  done  in  graphical  prognostic 
techniques  to  derive  possible  predictors. 

(c)  experimentation  with  the  use  of  absolute  vorticity  as  the  predictand 
and  recovering  the  height  field  by  relaxation  methods. 

(d)  experimentation  with  the  use  of  the  height  gradient  as  a  predictand 
rather  than  the  value  of  the  height  at  a  point. 

(e)  analysis  of  results  incorporating  operational  prognoses  rather  than 
perfect  prognoses. 

If  the  necessary  effort  could  be  applied,  implementable  results  could  be  made  available 
in  nine  to  twelve  months. 
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APPENDIX 


PREDICTION  EQUATIONS 

The  prediction  equations  derived  from  the  regression  analysis  have  the  form 
y  =  Ao  +  AA  +  a2x2  +  ...  -  Anxn,  (A-l) 

where  Y  is  the  predictand,  the  A’s  are  constant  coefficients  derived  (by  the  method  of 
least  squares)  from  the  developmental  sample,  and  the  X’s  are  the  predictors  selected 
by  the  screening  procedure. 

Each  set  of  prediction  equations  consists  of  six  equations:  the  three  predictands 
of  100-,  50-,  and  30-mb  height  change  for  24  and  48  hr.  Equations  from  experiments 
5  and  6  have  not  been  included. 

The  pair  of  numbers  that  is  associated  with  a  given  predictor  in  the  equations 
refers  to  the  grid  location  in  the  (K,  L)-grid  system  of  Figure  2.  The  symbols  and 
units  used  are  defined  in  Tables  I  and  II.  Note  that  while  the  error  statistics  in  the 
report  are  in  whole  feet,  the  predictands  in  these  equations  are  in  tens  of  feet. 
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Exp.  1  (north-south  orientation,  no  prognoses) 
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+0-2897  Z200(1>5)  -  °-2206H100-30(3>3) 
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z  (48)  =  -799.91  +  0.5666  AZ50(1,  3)  +  0.3956  AZ^l,  5)  -  0.0543  H100_3()(l,  3) 

-0.0777  Z50(5,  5)  +  0.1559  Z200(l,l)  -  0.3213  Z5Q(3,  3)  +  0.2550  Z30(5, 1) 
+0.3634Z50(1,5)  -  0.4277  Z5()0(l,  5)  -  0.4645^^(1,5) 

Z  (48)  =  -183.59  +  0.3662  AZ30(1,  3)  -  0.5136  H100_3()(3,  3)  +  0.2948  AZ^l,  5) 

+0.5383  AZ50(3,  3)  -  0.0774  Z5Q(5,  5)  +  0.0276  Z^l,  1)  -  0.3107  Z2Q()(3, 3) 

+0.7395  Z30(5,l)  +  0.1545  H100  30(l,  5)  -  0.5978  Z^,  1) 

-M91’100W3’3)  +  0.5759  H500_300(l,  5)  -  0.2361  Z3{)0(1,  5) 

+0.1536  Z  (1,  3)  +  4, 529, 800  Atj  (3,3) 
o  U  ou 
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Exp.  2  (north-south  orientation,  with  prognoses) 


A^100(24)  =  "3-3696  +  0.5266  AZ200(24)(3,  3)  +  0.2097  AZ^l,  3)  -  0.3916  ^200_100(3, 3) 

+°'1808H200-100(1’3)  +  °-2000H100-30(3>3)  -  °-1270H100-30(1’3) 

+0.2589  AZ50(3,  3)  +  0.0572  AZ200<24)(5f  3)  -  °-0175  H100_30(5,  5) 

+0.0541  AH100  30(1,  5)  +  2,060,500  At]5o(3,  3)  +  0.0462  AZ20Q(24)(1,  1) 

+0.0850  AZ500(24)(5,1)  +  0.10Q3  AZ^l,  1)  -  0.0500  AZ50Q(3, 3) 

+0-1096H500-300(1’5)  -  °*°696Z500(1’5)  "  °*0741Z500(3’3) 

+0.0452  Z5Q(1,  5) 

AZ  (24)  =  -479.07  +  0.2889  AZ  (1,3)  +  0.2413  AZonA(24)(3, 3)  +  0.2926  AZoa(3,  3) 

oU  oU  200  30 

-°-°141H100-30(3’3)  +  0-0602  AZ2oo(24)(5,  3)  +  0.0519  AH^  ^l,  5) 
+2,598,300  Ar?5()(3,  3)  +  0.1470  AZ3Q(1, 3)  +  0.0432  AZ20Q(24)(5,  5) 

+0.°954  Az5oo(24)(5»  !)  +  0.1354  H200_100(l,  3)  +  0.1198  Z^l,  1) 
-0.1916Z50(3,3)  +  0.1705  H500  300(3, 3)  +  0.1039  AZ^Jd,  1) 

+0.1197  Z5Q(1,  5)  +  0.0985  Z3q(5,1)  -  0.1263  Z50()(l,  5) 

+2,036,900  Ar?30(3, 3)  -  0.1131  5) 

AZ30(24)  =  “265*73  +  0.4563  AZ50(3,  3)  +  0.3279  AZ3q(1,  3)  -  0.2443  H1Q()  (3,3) 

+0.0294H2oo-ioo(3’3)  +  3,001,900Aij5()(3,3)  +  0.1338  AZ20Q(24)(3,  3) 

+0.1126  H100  30(l,  5)  +  2,780,400  Ar}3()(3, 3)  -  0.0440  Z5Q(5,  5) 

+O.2567  z30(5,  i)  +  0.0983  AZ5qo(24)(5,  1)  -  0.1447  Z5Q()(3,  3) 

-315‘2507J300(3*3)  "  °-1786Z50(5’1)  +  °-1060H200-100(1’3) 

+0.1286  AZ5oo(24)(1,1)  +  0.0958  Z^l,  1)  +  0.0737  AZ500(24)(5,  3) 

+0.1046  AZ100(5,  5)  +  0.0940  AZ500(5,  3) 
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AZioo(48)  =  -505.36  +  0.7019  AZ20Q(48)(3, 3)  -  0.6575  H200_10()(3,  3) 

+0.2482AZ50(1,5)  +  0.0771  Az200(48)(5,  3)  +  0.1680  AZ^l,  3) 

+0-3106H200-100(1-3)  +  °'0848H100-30(3'3)  '  0  2473  Z500<3'3) 
+0.1178Z300(1,3)  -  0.1617aZ5|)0(48)(3,3)  +  0.0850  Z3()(5, 1) 

+0-0876H200-100(1’5)  +  °'2307H500-300(1’5)  “  °'0938  Z500(1’ 5) 

+0.0997  AZ500(48)(5,  1)  +  0.1345  AZ500(24)(1,  1) 

AZ  (48)  =  -227.75  +  0.2911  AZ5()(1,  3)  +  0.2702  AZ20Q(48)(3,  3)  +  0.0086  H100_3()(3» 3) 

+0.2956  AZ30(1,  5)  +  0.1412  AZ2qo(48)(5,  3)  +  0.3114  AZ^l,  3) 

+0.2067  AZ&oo(48)(5,  1)  +  0.3145  AZ50Q(24)(1,  1)  -  0.0637  AZ200(1,  3) 
-0.4746  Z50Q(1,  5)  +  0.2482  Z300(l,  5)  +  0.2393  Z^l,  1)  -  0.4028  Z5Q(3,  3) 
+0.0962  Z50(5,l)  +  0.1575  Z50(l,  5)  +  0.2729  H500_300(3> 3) 

-0.1199  Z100(5,  5)  +  0.1717  Z100<5, 3)  +  0.1837  H^^l,  3) 

AZ  (48)  =  -699.10  +  0.4809  AZ  (1,3)  -  0.3604  H10()_30(3,  3)  +  0.1889  AZ2oq(48)(3,  3) 

30  o  U 

+0.1936  AH100  30(1,  5)  +  0.2750  AZ5()(3,  3)  +  0.2185  AZ50Q(48)(5,  1) 

-°-1520Z500(1’5)  +  °-3754H500-300(1>5)  “  0-0793  Z5()0(5,  5) 

+0.5851  Z30(5,l)  +  0.1228  AZ500(48)(5,  3)  +  0.1265  H200_10()(l,  3) 

+0.2281  AZ500(24)(1,1)  -  0.2195  H100_3Q(1,  1)  -  0.1548  Z5Q(5,  3) 

+0.1577  H100_30(l,  5)  -  0.1470  Z500(3, 3)  +  0.2462  AZ^l,  5) 

+0.2491  AZcn(5,  3)  -  0.2573  Z  (5,  1) 

50  oU 
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Exp.  3  (flow  orientation,  no  prognoses) 


AZ  (24)  =  -9.3100  +  0.2207  AZfl,  3)  +  0.3918  AZi3, 3)  -  0.0540  Z  (1,3) 
100  100  500  oU 


+0.1927  AZon(3,  3)  +  0.0943  Z3Q(5,  3)  +  0.1243  2^(3,3)  +  0.0502  Z[_nn(5,  5) 


30 


50' 


500' 


+2,785,800  Ar)c^n(3,  3)  +  0.2819  AZ^l,  3)  -  0.3645  nn(3, 3) 


'500 


100' 


+0.0965  Z  (5,  3)  +  0.0764  Zi5, 1)  +  0.2451  H, 


(3,3) 


'300'"’“'  - 500'“’  500-300 

+1,755,100ST500_300(3,3)  +  0.1703  AZ^l,  3)  -  0.0751  Z^l,  3) 


200' 


+0-1004H500-300(1’5) 


AZ5o(24)  =  73.1790  +  0.3935  AZ5()(1,  3)  +  0.4716  AZlQ()(3,  3)  +  0.2734  AZ3Q(1,  3) 

+0.2334  AH100  30(3, 3)  -  0.0138  Z3Q(1,  3)  +  0.0887  Ziq()(5,  1) 

+3,157,300  At?50(3, 3)  -  0.0665  Z5Q(3,  3) 

AZ30(24)  =  -282.83  +  0.5481  AZ5Q(3,  3)  +  0.3406  AZ3Q(1,  3)  -  0.2432  H100_3()(3,  3) 

+0.1504  H200_100(3, 3)  +  0.2431  Z30(5,l)  +  3,174,300  Ar^,  3) 

+2,494,700  At)3()(3,  3)  -  0.0085  Z2Q()(1,  5)  +  0.1460  AZ3Q(5,  3) 

+0.2317  AZ50(1,  3)  -  0.1609  Z50(5,l)  +  0.0606  H100_3()(l,  5) 

-0.0499  Z100(3,  3) 

AZioo(48)  =  284,40  +  °-2774AZ10o(1’3)  +  °-0352H200-100(1,5)  +  O-3235^1’5) 

-0.4972  Z10q(3,  3)  -  285,230^(3,  3)  +  0.2211  Z.^5,  3)  +  0.2131  2^(5,  3) 

+0.4454  AZ^l,  3)  -  0.1491  AH200.100(3, 3)  +  0.  5644  Z^S,  1) 

+0-0754Z500(5'5)  +  0'2025AZ200(1'3)  +  °'2219  H500-300<3' 3) 

+0.0663  Z200(l,  5)  -  0.6595  Z300(5,l)  +  0.3910  Z20Q(5, 1)  -  0.1908  Zl0Q(5, 3) 
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AZ50(48)  =  -246.20  +  0.3621  AZ^l,  3)  -  0.1228  H100_3()(l,  3)  +  0.4462  AZ^l,  3) 

+0.1763  AZ  (1,5)  -  0.3984  Z50(3,  3)  +  0.1326  Z1Q()(5,  1) 

+0.1949  Z30(5,  3)  -  0.2248  Z500(l,  3)  +  0.0777  Z^,  3) 

+0.2040  AZ30(3, 3)  +  0.1630  Z5Q(1,  3)  +  0.2064  AZ5Q(1,  5) 

+0.1946  AZ500(1,  3)  +  0.0542  Z300(3,  3) 

AZ  (48)  =  -142.22  +  0.5430  AZ^l,  3)  -  0.5100 ^^  ^(3,3)  +  0.6605  AZ5Q(3, 3) 

+0.3957  Hi00  3()(5,  3)  +  0.2485  AZ^l,  5)  -  0.3572  Z^,  3) 

+0.2128  Z10q(5,1)  +  5,077,600At?5o(3,3)  -  0.0824  Z^l,  5) 
-877.57°,500(3,3)  +  0.1108  H200  100(l,l)  -  0.3335  3) 

+0.3071  AZg0(l,  3)  +  0.2880  H200  100(l,  3)  -  0.3151  1) 

-0.5212  z2oo<3’  3)  +  °-2267H500-300(1'3)  +  °'4773  Z300(3' 3) 

+0.4039  Z  (5,3)  -  0.2766  Z5q(5,  3) 
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Exp.  4  (flow  orientation,  with  prognoses) 


AZ  (24)  =  -79.455  +  0.5704  AZ200(24)(3, 3)  +  0.2455  AZ^l,  3) 

+0.0747  AZ200(24)(S,  3)  +  827,690ST,,00_100(3,  3)  -  0.0015 H100_30<1,  5) 
+0.1719AZ50(3,3)  +  0.1003  AZ100(1, 3)  +  0.0877  AH100  3(|(3, 3) 

-°-3520H200-100(3'3)  +  °-°837H100-30(3'3)  +  °-1836H200-100(1’3) 
-0.0696  Hi00.3o(1.  3)  *  0.1542  Z500(3'  3)  *  °'0484  Z300<5’ 3) 

-0.1144  AZ500(24)(3, 3)  +  0.0546  Z3C<5, 3)  +  0.0442AZ20(|(24)(1,  1) 

AZ  (24)  =  151.55  +  0.3356  AZ  (1,3)  +  0.2117AZ  (24)(3,3)  +  0.2597  AZ  (3,3) 

50  50  i^UU  ou 

-0.0621H100  30(3,3)  +  0.0701  AZ2qo(24)(5,  3)  +  0.2619  AZ^l,  3) 
+3,229,900  A775o(3,  3)  +  0.0555  AZ2oq(24)  (5,  5)  +  0.1261  AH^  ^l,  3) 

AZ30(24)  =  -32.697  +  0.5201  AZ5Q(3,  3)  +  0.4367  AZ3()(1,  3)  -  0.2047  H100_3()(3,  3) 

+0.0114  H200  100(3,  3)  +  0.0713  Z30(5,l)  +  3,128,200 At)5()(3,  3) 
+2,865,200Arj3()(3,3)  -  0.0411  Z^l,  5)  +  0.1082  AZ.^24) (3,  3) 

+°.1603H200-100(1’3)  "  445’110W3’3)  -  °-0978Z500(3>3) 
+0-0751H200-100(5’5)  +  0-1265  AZ50(5,  5) 


AZioo(48)  =  ?6-952  +  0-7417  AZ2oo(48)(3,  3)  -  0.6964  H200_10()(3,  3) 

+0.3551  AZ50(1,  3)  +  0.0248  AZ20Q(48)(5,  3)  +  0.2477  AZ5()(1,  5) 
+0.0830  H100  30(5,  3)  -  0.3246  Z500(3,  3)  +  0.0213  Z^l,  3) 

-0.1980  AZ500(48)(3,  3)  +  0.1172  Z5Q(3,  3)  +  0.0616  AZ20()(48)  (5,  5) 
+0.0833  AZ200(48)(5,1)  +  0.1082  Z500(5, 1)  +  0.1199  H200  10()(l,  3) 
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AZrQ(48)  =  334.61  +  0.5992  AZ50(1,  3)  +  0.3150  AZ20()(48)(3,  3) 

-0.1658  H100_30(3,  3)  +  0.3219  AH100_30(1,  3)  +  0.09S2  AZ200(48)(5,  5) 
t0.1031H1(|0_30(5,3)  +  0.2543  AZ3q(1,  5)  +  0.2423  AZ3(|(3,  3) 

+0.1359  AZ200(48)(5, 1)  +  3,237,900  Al^fS,  3)  +  0.0703  AZ20()(48)(1,  1) 
+0.1521  AH200  100(1, 3)  -  0.0651  V(3, 3)  -  0.1230  Z500<1,  5) 

+0-1586H500-300(1’5)  +  °'1225Z200(5’1)  ‘  0  0945  Z50(3':i) 

AZ„  (48)  =  543.78  +  0.5533  AZ3()(1,  3)  -  0.3323  H100_3()(3, 3)  +  0.2203  AZ20()(48)(3,  3) 

3  0 

+0.5221  AZ5Q(3,  3)  -  0.1457  Z50Q(1,  5)  +  0.1606  Z3()(5,  1) 

+0.1849  AH100_30(1,  5)  -  0.0691  Z500(5,  5)  +  0.1023  AZ20()(48)(5(  3) 

-0.2079  AZ900(3, 3)  +  3,834,700  At?5()(3,  3)  -  0.1483  AZ20()(1,  3) 

+0.3249  AZ50(1,  3)  +  0.0900  AZ50()(48)(1,  1)  +  0.1923  AZ5()(1,  5) 

-0.0956  Z500(3,  3)  -  744,240t)2()0(3,3)  +  0.1706  H200_10()(l,  3) 

-°-2627H500-300{1’1)  '  °-1172H100-30(1*1) 
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